Electron energy spectrum of the spin-liquid state in a frustrated Hubbard model 
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Non-local correlation effects in the half-filled Hubbard model on an isotropic triangular lattice are 
studied within a spin polarized extension of the dual fermion approach. A competition between the 
antiferromagnetic non-coUinear and the spin liquid states is strongly enhanced by an incorporation 
of a k-dependent self-energy beyond the local dynamical mean-field theory. The dual fermion correc- 
tions drastically decrease the energy of a spin liquid state while leaving the non-coUinear magnetic 
states almost non-affected. This makes the spin liquid to become a preferable state in a certain 
interval of interaction strength of an order of the magnitude of a bandwidth. The spectral function 
of the spin-liquid Mott insulator is determined by a formation of local singlets which results in the 
energy gap of about twice larger than that of the 120 degrees antiferromagnetic Neel state. 

PACS numbers: 31.15.V, 71.10.Fd, 75.10.Jm 



I. INTRODUCTION 

In some special cases, quantum fluctuations can de- 
stroy an ordered state even at zero temperature. The 
study of such "quantum phase transitions" due to zero- 
point fluctuations is a hot subject in modern condensed 
matter physicJ^'. These fluctuations are especially im- 
portant in frustrated systems. A prototype of classical 
frustrated spin system is the Ising model on a triangular 
lattice with nearest-neighbor antiferromagnetic interac- 
tion. The ground state of this model has a macroscopi- 
caUy large degeneracy and the addition of quantum spin- 
flip terms results into a formation of spin-liquid state 
known as "resonating valence bonds" (RVB^P^l There 
are strong experimental evidences that the RVB is the 
ground state of some Mott insulators: organic salt k- 
(ET)2Cu2(CN):^and spinel oxide on the frustrated hy- 
perkagome lattice Na4lr30^. 

Theoretical description of the materials with the RVB 
Mott insulating state requires an accurate choice of pa- 
rameters of the model. The Hubbard model on a trian- 
gular lattice with nearest-neighbor hopping is the sim- 
plest choice, though its validity can be questioned. First, 
recent ab-initio calculation^^ suggest that the family 
of K — {CN) Mott insulating organic salts has a large 
anisotropy of approximately 20% difference between hop- 
ping parameters in different directions. Moreover, it 
is not clear which phase of the Hubbard model on a 
isotropic triangular lattice is a ground state. The RVB 
phase and 120 degrees antiferromagnetic Neel states are 
expected to possess very close energie^i^llil. Any of 
these phases can be stabilized in the experiment by tun- 
ing next-neighbor hopping, electron-phonon interaction, 
anisotropy etc. This makes the results for the ground 
state obtained within the Hubbard model on an isotropic 
triangular lattice not straightforwardly comparable with 
an experiment. 

In the limit of a strong on-site repulsion U this model 
is equivalent to the Heisenberg model with antiferromag- 



netic nearest-neighbor exchange interaction J — 2t^ /U 
where t is the hopping parameter. The classical Heisen- 
berg model predicts a noncollinear magnetic structure 
with 120° angles between neighboring spinJl^l, Prop- 
erties of quantum Heisenberg model are likely to be 
similar and several theoretical and numerical QMC and 
DMRG results predict the Neel 120° ordering of ground 
statdJ^tll^. It was suggested that a disordered state of 
quantum Heisenberg antiferromagnet may behave as a 
spin liquicP^. 

For moderate values of U , the situation is less clear. 
The results o btain ed by a path integral renormalization 
group (PIRG/^SE^ and the variational cluster approach 
(VCA)^^ suggest the non-magnetic insulating state to 
prevail at zero temperatures in certain interval of pa- 
rameter U/t and the Neel ordered state to exist at higher 
values of U. 

To study the Hubbard model for the strongly corre- 
lated lattice systems, the dynamical mean-field theory 
(DMFT)-^Q''22]is frequently used. This ap proach has been 
employed to the case of triangular lattic^^HSll^ however, 
the localDMFT scheme is not sufficient to study the RVB 
physics based on the formation of sing lets on bondtPEHH 
which involve non-locality, at least, at the level of pairs 
of sites. Recently, non-local methods beyond the DMFT 
such as the cluster extension^ of DMFT dynamical clus- 
ter approximation (DCAjf^IMl and dual fermions (DF) 
theorji^^ have been applied to the problenPiHlIl How- 
ever, the competition of the two phases has not been con- 
sidered there and their electron spectral function phases 
has not been discussed. 

In this paper a spin-polarized extension of the DF tech- 
nique is implemented for the isotropic triangular lattice. 
Contrary to the local DMFT, the long-ranged electron 
correlations are taken into account within a regular dia- 
grammatic perturbation scheme starting from the DMFT 
as zero-order approximation^^. It is shown that these 
corrections change essentially the phase diagram of a tri- 
angular lattice decreasing the energy of the spin-liquid 
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state in comparison with the DMFT results. For the 
120° Neel State the corrections are much smaller so that 
the spin liquid state has a lower total energy than the 
Neel-ordered phase in a certain range of parameters. We 
calculate the electron energy spectrum by means of the 
DP scheme utilizing the continuous time quantum Monte 
Carlo solveJsS at low temperatures. 



Rm is a radius vector of the supercell, points to a 
site inside it. K is defined inside a Brillouin zone of 
a supercell, k is its basis vector. In this notation, the 
Hamiltonian of the Hubbard model can be written as 
follows: 



H ^ 



4K^(£kk'(K) - Ai)ck'Ko 



[Ufii^hii, (3) 



Kkk'cT 



II. METHOD 

In the spirit of a mean-field theory the whole lattice 
problem is mapped to the single impurity, in this case 
to a single-impurity Anderson model subject to a self- 
consistent condition. Such a dynamical mean-field the- 
ory (DMFT) restricts an investigation to studying only 
local correlation effects. For frustrated quantum systems 
this approximation is obviously not applicable. In or- 
der to take into account non-local corrections two ap- 
proaches are implemented: (i) the long-ranged correla- 
tions for the anti-ferromagnetic ordering are included by 
introducing the spin polarization to the effective bath, 
(ii) the short-ranged correlations are incorporated within 
the dual fermion technique - it is based on a transfor- 
mation to new dual fermionic variables, which make the 
diagra mmati c expansion over non-linearity much more 
effective™. 

We start with a single-site DMFT-like problem. For a 
given impurity with a yet undefined dynamical hybridiza- 
tion function A^^o-, which describes the interaction with 
a fermionic bath of surrounding lattice, the effective ac- 
tion for single-impurity Anderson model is written in the 
following way: 

( A^,^ - A* - «^)cLctJT + / n^T-nirdT, (1) 
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where index a ^ labels the spin projections, U is an 
a on-site Coulomb interaction, /i is a chemical poten- 
tial, and /3 is an inverse temperature. The solution of 
this problem can be obtained numerically by the CT- 
QMC solvei'^. Along with the one-particle Green's func- 
tion = —i{Ci^<jcl,a) one can obtain the higher or- 
der correlators, such as the two-particle Green's function 
X1234 = (Tco-i(Ti)ccr2(T-2)43('^3)44('^4)) and the corre- 
sponding vertex function 

7^4 = .9n'.922^ [Xl'2'3'4' - (.9l'3'g2'4' " 5l'4'52'3' )] 5^334^4 

We now need to establish the connection between the 
single-site model and original lattice problem. The 
whole lattice is divided into supercells, combined of three 
atoms. The Fourier transforms over both the supercell 
and cluster spatial indices are utilized, which corresponds 
to the DCA-like schem^ 

G(r = Ra/ + n, r' = Rm' + rj) 

Gkk'(K)e'K(^'"-^^^')e^('''"'~'''""j), (2) 



where index i labels a of site of the lattice. For the super- 
cell approach the dispersion law is given by £kk'(K) — 
e(K -|- k)5kk' with the spectrum of a triangular lattice: 



r 2tt 

e(K = aki + /^ka) = -2t cos( — (-a + (3)) + 

> 27r , , , , 27r , 

-I- cos( — (a 2/3)) + cos( — (-2a - /3)) 



(4) 



where ki , ka are non-coUinear basis vectors of the su- 
percell in the reciprocal space: ki = |(kx -I- ky), ka = 
|(— kx -I- 2ky),|ki| = jkal, where \<i.xi^y are the basis 
vectors of the reciprocal triangular lattice. 

After adding and subtracting an arbitrary hybridiza- 
tion function into the action for the Hubbard model, one 
can obtain: 

5[c, c*] = E c^ck.(A.kk' -ekk'(K))cKkv 

(5) 

It should be mentioned that the k dependence of A^^ 
arises from the spatial Fourier transform over internal 
sites of the supercell. In case of all sites being equivalent 

Ai^kk' — A^Jkk'- 

The further steps of the dual fermion scheme for a 
supercell are similar to the cluster DF-approaclP^. By 
utilizing Hubbard-Stratonovich transformation to new 
fermionic variables and integrating out the original de- 
grees of freedom one can obtain a new fermionic problem, 
which is even more complex that the original one, but it 
allows to be efficiently treated by a diagrammatic expan- 
sion by a right definition of a yet arbitrary hybridization 
function ts.^^. It is chosen in a way for the trivial case of 
no diagrammatic expansion of the self energy of the dual 
variables to be equivalent to the DMFT. We use a first 
non-zero term in the diagrammatic expansion of a dual 
self- energy: 



yd _ 
^11' — 



1 



"lji+W3 7i234 



G'44/(R)G'2'2( 



-R)Gf,3(R)7if^,2'i' 



(6) 



K 
kk' 



The important long ranged corrections to the single-site 
DMFT scheme are introduced in formula ([6|. It should 
be noted that the dual fermion Green's function con- 
sists only of non-local part while the vertex function 
is local and obtained within a single-impurity Anderson 
model^^ 
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In order to introduce non-coUinear spin-dependent DF- 
scheme one can rotate spinors for different sites of a su- 
percell. The local Green's function for a given site is 
obtained as follows: 

TT /IT - TT yTT 

g^,^U,{-,±-)g^U:{-,±-), (7) 

where [/ is a spin 2x2 rotation matrix. 

The total energy (H — iiN) for different states has been 
calculated in the following fornP^ 

{H) = {Ho + V-^iN), (8) 
{Ho}=-tJ2Goj{T = 0+), (9) 

{H, + 2V-f^N)^(^^^) (10) 

This procedure requires the Fourier transform of the 
Green's function to the imaginary time domain. Since 
the CT-QMC solver provides a bare Green's function on 
finite number of Matsubara frequencies, this could be a 
serious source of numerical errors. The interaction ex- 
pansion solveil^Sl produces low noise at the high Matsub- 
ara frequencies and allows to fit an asymptotic tail of the 
Green's Function. 



III. RESULTS 

In order to calculate the energetics of a triangular lat- 
tice, 2-3 different independent runs consisting of 20-30 
DMFT iterations and following 15-25 DF iterations were 
done with the lattice size of 32 x 32. The interaction 
expansion QMC solver'^ was used for which the QMC 
error becomes sufhcient at high values of U. The ob- 
tained data was averaged, and the corresponding statis- 
tical error is plotted on the graphs. The spin-polarized 
phase in the insulating regime hence being constructed 
by spin rotation represents the Neel state. The unpolar- 
ized DMFT insulating phase has only local correlations 
and does not possess any magnetic ordering. It will be 
shown that the non-polarized insulating DF phase is con- 
structed by singlet corrections to a local self-energy and 
therefore perturbatively describes the spin-liquid state. 

A. Total energy calculations 

The results for a total energy per site {H + iiN) = 
{H) + U/2 and averaged square of magnetic moment (m?) 
as a function of U/t in the region of 8 < U/t < 14 for 
various states at half-filling, /3t = 10 are presented at 
Figures 1 and 2 respectively. The separate curves for the 
phase transition from a metallic to a spin liquid state and 
from a metallic to a 120° Neel state were obtained. The 
phase transition to a spin liquid state occurs at U/t — 
9.6±0.2, the hysteresis a.t U/t = 8.25— 9±0. 2 represents a 




U/t 



Figure 1. (Color online) The dependence of the total energy 
per site at /3t = 10 on the parameter U. The separate curves 
for a phase transition from metallic to a spin liquid and from 
metallic to a Neel state were obtained. The spin liquid state 
possesses lower energy in the interval 9.5 < U/t < 13. At 
higher values of U the energies of these insulating states be- 
come close. The phase transition from metallic to the spin 
liquid phase occurs at U/t ~ 9.6 ± 0.2, the value of the criti- 
cal U/t for the Mott transition from metallic to the 120° Neel 
phase lies in the interval 8.25 — 9.4±0.2. The dashed line rep- 
resents the monotonic interpolation of the spin liquid energy 
curve. 
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Figure 2. (Color online) The dependence of the averaged 
square of magnetic moment per site at j3t = 10 on the 
parameter U. With the increasing of U the formation of 
magnetic state is detected first for the 120° Neel state at 
U/t = 8.25 - 9.4 ± 0.2. The hysteresis in the Neel curve rep- 
resents a first-order phase transition from metallic to a Neel 
state. The formation of the magnetic moments for the spin 
liquid state occurs at the value of critical [/ of a Mott phase 
transition (7 = 9.6 ± 0.2. 



first order transition from the metallic to the Neel ordered 
state. 

At first, it can be seen from the energy graph that 
there is an essential difference between the DMFT and 
DF results for the spin liquid case. Non-local correc- 
tions favour the unordered insulating state. Without DF- 
corrections this phase never corresponds to the ground 
state in the whole studied region and remains metallic up 
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to U/t < 11.2 ±0.2. This is in agreement with the previ- 
ous DMFT calculation^^. Thus the first non-vanishing 
diagram in the DF approach lowers the total energy of 
the unordered insulating state to make it favourable in 
the region of 9.5 < U/t < 13. The dual corrections in the 
Neel phase are small and do not change the behaviour of 
this phase. At higher values of U/t the curves of total 
energies of Neel and spin liquid state become close and 
cannot be distinguished within the statistical error of the 
method. Nevertheless, since the dual corrections vanish 
at high Of22l one would expect that the total energies of 
DMFT and DF curve will become equal, with the Neel 
ordering curve to dominate in this region. 

A formation of local magnetic moments with the in- 
creasing of U/t is captured in the Figure 2 {f3t — 10). It 
can be seen that in the insulating region both spin liq- 
uid and 120° Neel states possess well-defined magnetic 
moments thus differing by a type of ordering. The dual 
correction for the unordered state significantly enhances 
a formation of magnetism in the region of Mott transi- 
tion, thus emphasizing the role of singlet correlations. As 
for the Neel phase these corrections are small and mag- 
netic ordering is formed even at lower value of U, which 
means that a static spin ordering is essentially captured 
on the DMFT level. 



B. The density of states and the nonlocal 
self-energy 

The densities of states N{E) for the strongly corre- 
lated metallic state are shown in the Figure 3. One can 
clearly see a three-peak structure typical for the systems 
at the metal side of Mott transition, namely, two Hub- 
bard bands and central "Kondo" quasiparticle peak^. 
This general structure is very similar both in the local 
DMFT and the non-local DF approximations. There is 
an essential difference between two schemes related to 
the shape of the Kondo peak. The density of states at 
the Fermi level in the local approximation is fixed by the 
sum rule, namely, N{0) should be equal to its value for 
noninteracting electron^^D, jn the DF approach, this is 
no more the case and one can see that the non-local cor- 
rections suppress the height of the Kondo peak at larger 
U. 

The densities of states for the Neel and spin liquid 
states are shown in the Figure 4. One can see that the 
spin liquid solution has a wider gap than the Neel state. 
This agrees with the RVB-picture as a liquid formed by 
local singlets^. Indeed, local singlet must have a spec- 
trum inherited from a spectrum of a two-site cluster. To 
emphasize this feature, we show a spectrum of the iso- 
lated 2-atom cluster with the same value of t with black 
ticks in the Fig. 4. The grey ticks represent value zLU/2. 
As one can see from the inset of the Fig. 4, the non- 
local corrections to nearest-neighbor self-energy in the 
spin-liquid case are much larger than in the Neel state 
and dominate the low-frequency spectrum. This reflects 
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Figure 3. (Color online) The density of states N{E) of a 
metallic phase at /3t = 10 and U/t = 8.8, 9.2, 9.4. The upper 
part of the graph corresponds to a dual fermion solution, the 
lower part is obtained by the DMFT. The height of the central 
"Kondo" peak calculated by a DF technique is not fixed. 
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Figure 4. (Color online) The density of states N{E) of the 
insulator phase at /3f = 15 and U/t = 11.2 for Neel and spin 
liquid states. The difference between insulating gaps is seen. 
The spectrum of the 2-atom singlet with the same values of 
U and t is plotted by black ticks. The grey ticks represent 
the value ikU/2. The inset shows the real part of the nearest- 
neighbor self-energy over Matsubara frequencies in spin liquid 
regime (black). The same value obtained for Neel state (red) 
is negligible, while non-local self-energy function in DMFT is 
by definition equal to zero. 



an importance of nonlocal dynamical fluctuations for the 
spin liquid state. On the contrary, for the Neel phase a 
most of the nonlocality is reflected in the static ordering 
and is therefore taken into account already at the DMFT 
level. This corresponds to the results of the Fig. 2. 



IV. DISCUSSION 

Some approximations have been made in our consid- 
eration, and they should be discussed in the context of 
the phase diagram of the model studied. The main of 
them are a finite-temperature consideration and taking 
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into account just the leading dual correction. 

The dual fermion corrections are, in general, sufficient 
in the intermediate values of U /t at the order of a mag- 
nitude of the bandwidth. The picture obtained by in- 
cluding first non-zero diagram corresponds to an idea of 
a spin-liquid as a state which can be observed in certain 
region of parameter U /t with metallic and Neel ordered 
states to occur outside of this rang e. This corresponds 
to what is provided in the Iiteratur^l2lini36l_ p^^ ^j^g Yig. 
2 it was shown that an inclusion of singlet corrections 
sufficiently enhances the magnetism and makes the for- 
mation of magnetic state to occur at lower values of U . 
Our calculations have been done using a leading diagram 
only and we expect higher order diagrams to shift the 
value of Uci even lower, so that a spin liquid state should 
have a monotonically increasing total energy curve - that 
is illustrated by a dashed line in the Fig. 1. We expect 
a spin liquid ground state to be extended to the region 
beyond current Uci = 9.6 and a Neel ordered state to 
appear after the spin liquid. To have more quantitative 
description of a phase diagram one has to include a se- 
ries of diagrams, such as ladder, which may lead to a 
formation of bounded sing let stated 

The upper critical value Uc2 > 14 looks overestimated 
due to a finite temperature. Indeed, thermal fluctua- 
tions in general destroy an ordering. The most clear case 
is a half-filled Hubbard lattice, where an ordered ground 
state is realised at any (even infinitesimal) C/, due to so- 
called perfect nestling. However, DMFT calculations at 
relatively small temperature = 20 show a phase dia- 
gram with the ordering destroyed by local fluctuations of 
magnetic momenta up to J7 w Ipflj. Speaking about the 



energy dependencies presented in the Fig. 1, we expect 
that lowering the temperature should shift the energy of 
the Neel state down. 



V. CONCLUSIONS 

In summary, the Hubbard model on a triangular lat- 
tice is investigated within a framework of spin-polarized 
dual fermion approach. The ordering of 120° Neel phase 
is described by an introduction of a non-collinear spin- 
polarization via a spinor rotation to the effective bath. 
The DMFT approximation works well for this type of 
ordering since non-local correlations are relatively small. 
On the contrary the solution with a non spin-polarized 
bath which describes the spin-liquid state has strong non- 
local corrections. This leads to a significant change in the 
phase diagram at low temperatures - the DF approach 
produces an insulating state with the lower total energy 
compared to that of a 120° state in a certain region of 
parameter U /t. Note that the simple DMFT solution 
gives an artificial Fermi liquid at these parameters. The 
densities of states for a spin-liquid phase has much larger 
insulating gap than for a non-collinear Neel state. Due 
to sufficient non-local corrections the dual fermion ap- 
proach can capture the essential physics of a spin-liquid 
phase. 
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